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A duality argument for the recently discovered chiral-nematic phase of the XY model in a tri- 
angular lattice is presented. We show that a new Ising variable naturally emerges in mapping the 
antiferromagnetic J\ — J2 classical XY spin Hamiltonian onto an appropriate Villain model on a 
triangular lattice. The new variable is the chirality degree of freedom, which exists in addition to 
the usual vortex variables, in the dual picture. Elementary excitations and the associated phase 
transition of the Ising degrees of freedom are discussed in some detail. 
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I. INTRODUCTION 

A description of the excitations of a classical ferromag- 
net in a two-dimensional lattice is based on two entities: 
spin waves, which are small fluctuations of the spin ori- 
entation from their ground-state direction, and vortices, 
which are topological structures of spins carrying a non- 
zero winding number. Among the two, the latter exci- 
tation is more effective at destroying the coherence of 
spins and ultimately drives the phase transition between 
a quasi-long-range ordered (QLRO) phase and a disor- 
dered, paramagnetic phase 1 . A convenient mathematical 
description of the phase transition in terms of the vor- 
tex variable is afforded by the mapping first discussed by 
Villain (Villain mapping) 2 . 

There are instances, however, where vortices do not 
exhaust all the relevant excitations for the phase transi- 
tion. Another variable, namely, chirality, often appears 
in models of classical spins with frustration 3 . With its 
symmetry being Ising-like, the chirality fluctuation can 
drive a second phase transition, apart from the one driven 
by vortex proliferation, of fsing universality class. A gen- 
eralization of the Villain mapping for spin models with 
frustration was given by Villain 3 , where it was demon- 
strated that the chirality emerges as a new, independent 
excitation mode of the model. 

As a most recent example of a chirality-driven phase 
transition, it was shown in the context of a generalized 
antiferromagnetic XY spin model on a triangular lattice 
that there can be a chirality-driven phase transition tak- 
ing place at temperatures well above the magnetic tran- 
sition temperature 4 . The situation stands in contrast 
with the case of standard antiferromagnetic XY model 
on the same lattice, which only reveals a tiny separation 
between the chirality order (T x ) and the magnetic order 
(Tkt) temperatures 5 . 

The aim of this paper is to place the observations made 
in Ref. 4 in the context of the duality mapping of Villain 
by identifying the proper Ising variables associated with 
the chirality. While Villain's original paper 3 considered 
a square lattice, we focus here on a triangular lattice, as 
was the case for the numerical work of Ref. 4. A more 
elaborate discussion of Villain's original model appeared 



in the work of Lee and Grinstein 6 , but their discussion 
was for a square lattice, and the issue of the chirality did 
not arise in their work. 



II. DUALITY MAPPING 

The model studied in Ref. 4 reads 

H = Ji cos((pij) + J 2 ^ cos(2(pij), (1) 

{V) (ij) 

where ipij is the angle difference ipi — ipj between near- 
est neighbors (ij). Focusing on a pair- wise interac- 
tion J 1 cos[(pij\ + J 2 cos[2</?jj], the minimum energy angle 
for a given (ij) bond is two-fold degenerate at ir ± A, 
A = cos _1 (Ji/4J 2 ), provided 4J 2 > J\. For 4J 2 < Ji, 
a single minimum-energy angle is achieved, and the du- 
ality mapping proceeds in the same manner as in the 
ferromagnetic case. 

An appropriate Villain model for the interaction hav- 
ing a double minima is given by the pairwise probability 
weight 

P[ ( Pij]= exp (~Y [¥ij- ,7T -At lj + 2iTn ij } 2 ^ . (2) 

We introduce two integer link fields: riy, which run from 
—00 to 00, and tij, which takes on two values t{j = ±1 
corresponding to two equivalent minima. Both fields 
change sign under the interchange of the lattice indices, 
iiji = —riij, and tji = —tij. K is related to the inverse 
temperature 1/T. Going through the standard Villain 
mapping gives an alternative expression containing only 
the first power of ipij : 

p [<Pij] = exp [~^K l io +ihj(<Pij-n-M ij ) \ , (3) 

' ij 1 1 ij 

where is another integer field running from —00 to 00. 
The partition function Z is obtained as the product of 
P[(pij\s over all nearest-neighbor links: 
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The exponential factor in the above equation is evalu- 
ated for each configuration {hj,Uj} and is summed over 
all possible configurations to give Z. Integration over the 
angle fa can be done first, which imposes current conser- 
vation at each site. In turn, the constraint can be solved 
by introducing a dual integer field hi, which is connected 
to kj by 
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hi — hi + Ei, 
hi+E 3 — hi, 

hl-E 2 +E 3 — hi + E 3 , 
hl+E 1 -E 2 +E 3 — hi-E 2 + E 3 , 
hl+E x -E 2 ~ hi + E x -E 2 + E 3 , 
hl+E! — hi + E ± -E 2 - 



(5) 



See Fig. 1 for the definitions of various labels. 
Now, we have the partition function 



Z = 
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where A is given as the sum over the dual links (IJ). 
The spin-wave contribution to Z has been dropped. On 
re-organizing the expression on the far right of A for each 
hi, we get an equivalent expression 



5^(fcj-/lj)(7r + Aty) 

(IJ) 



J2hi(3ir + Aa). (7) 



Here, cj is the sum of the tjj's for each triangle cen- 
tered at / and going in the counter-clockwise sense. The 
allowed values of c/ are cj = ±3, ±1. The partition func- 
tion is reduced to 



Z= ex P (-^£(/ l/ -/ lJ ) 2 -^/ l ^+Ac / )] 

{hi, Uj} \ (U) I J 

(8) 

At this point, we invoke the Poisson sum formula to re- 
write the partition function as 
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Here, mi is another integer field from —00 to 00 defined 
at the dual sites /, which will play the role of the vorticity. 
We will integrate out Oi to obtain the partition function 
solely in terms of the vorticity mi and chirality cj : 



Z ~ ^2 ex P ( -2k 2 K MiGijMj j . 

{mi,tij} \ U ) 



(10) 



In the sum, I and J independently run over all dual- 
lattice sites. The real-space Green's function, Gu, is 
given by 
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(27r) 4 — 2 cos(fc x a) — 2 cos(k y a) ' 



(11) 



It is convenient to use the regularized Green's function 
Gu - G , so 
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Z = 



exp (-2n 2 KG (j2 Ml ) ) 
x exp (^-2it 2 K J^MiG'ijM,^ . 



(12) 



Since Go ~ (l/2n) \n(L/a), the allowed configurations 
are those with total vorticity J2i Mi = 0. With the help 
of the result at large distances 
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the partition function may be rewritten as 

/ 7T 2 

Z= expl K - — Kj2Mj+irKj2MiMj 

{mi,tij} I I^J 



In 
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(14) 



We have used the relation (£j Mi) 2 = = J2i M l + 
zC/^j MiMj. Often the fugacity term is defined by 
y = cxp(— n 2 K/2). This term controls the number of 
(total) vortices in the system. Equation (14) is the de- 
sired action expressed solely in terms of the vorticity m; 
and the chirality cj (or tij). 
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TABLE I: Local chirality cj and local vorticity mj for which 
the net vorticity Mi is least costly. 



In Rcf. 4, the large J2/J1 region was shown to have 
a chirality ordering transition taking place well above 
the magnetic transition. In the J2/J1 3> 1 region, A = 
cos~ 1 (Ji/4J2) becomes close to it/2, and one can write 
A = 7r(l — (5)/2, where <5 is given by sin(7T(5/2) = J1/4J2. 
The fugacity consideration requires that M| be as small 
as possible at low temperatures, and in Table I, the small- 
est net vorticity Mj for each chirality value cj is shown. 
Furthermore, Mf for c/ = ±3 is less than Mj when 

a = ±1, 



Mj(ci = ±3) - Mf( Cl = ±1) « -(5/2 < 0. (15) 

The global minimum of ^ 7 M| is obtained if each tri- 
angle carries the chirality a = ±3. In practice, this is 
possible by arranging all the up triangles to carry c/ = 3 
and all the down triangles to carry cj = —3, or vice versa. 
This is the chirality-ordered phase at low temperature. 

The cheapest excitation is the one with the least in- 
crease in Mf. Such an excitation is achieved if for a 
given triangle, two of the t^'s reverse their directions, 
and cj = 3 becomes cj = -1 or q = —3 becomes 
cj = 1. As one can see from Table I, the corresponding 
change in the net vorticity is from M/ = —1/4 + 35/4 to 
Mi = -1/4 -(5/4, or AM/ = -5. Similarly cj = -3 -> 1 
incurs the change AM/ = <5. Such excitations are called 
incommensurate vortices 7 to set them apart from integer 
changes in the vorticity, also called commensurate vor- 
tices. At a temperature 7j that scales with (5, we thus 
expect an Ising transition caused by the proliferation of 
incommensurate vortices. Imagine we draw a line seg- 
ment from I and J that intersects the (ij) bond with its 
tij reversed. Then as one can see from Fig. 2, the cre- 
ation of incommensurate vortices is achieved if the line 
segments close onto themselves, forming closed loops. In 
the dual hexagonal lattice, the smallest such loop would 
be a single hexagon. At the Ising transition, the size of 
the closed loop diverges. In the numerical study of Eq. 
(1), an Ising transition with T/ ~ 6 was identified for the 
large J2/J1 region. 

A flip of a single tij in a triangle, on the other hand, 



changes c/ = 3 (-3) to c/ = 1 (-1), and the net vortic- 
ity change AM/ = ±(1/2 - (5/2). Close to 5 = 0, this 
is creating an isolated half-integer vortex with a rather 
high energy cost. To avoid the creation of a half- integer 
vortex, all triangles must have Ac/ = ±4, rather than 
Ac/ = ±2. For an open line segment, if both ends of the 
line terminate at the A or the B sublattice sites of the 
hexagonal lattice, the endpoints are associated with half- 
integer vortices of the same sign. If one end terminates 
on a sublattice different from that of the other end, then 
one has created a pair of half-integer vortices of opposite 
charges. It is then the proliferation of closed loops in 
the dual hexagonal space that corresponds to the Ising 
transition of the chirality. At exactly 6 — (J\ = 0), 
the energy cost for the flip of chirality by four units be- 
comes zero, so already at zero temperature, one has a 
disordered chirality phase. At small 8, a proliferation of 
the cluster of triangles with Ac/ = ±4 or a proliferation 
of closed loops in the dual hexagonal space occurs at a 
small finite temperature Tj. On the other hand, once the 
length of the line becomes infinitely long, the half-integer 
vortices becomes free, and the half-integer vortex pairs 
unbind. This is the KT transition of half-integer vortices 
at a higher temperature Tkt- 



III. SUMMARY 

To summarize, we have presented a duality argument 
for the antiferromagnetic J\ — J2 XY spin model on a 
triangular lattice. In carrying out the Villain mapping 
of the original spin model, an additional Ising degree of 
freedom emerges naturally. This new degree of freedom 
corresponds to the spin chirality, and its excitation is re- 
sponsible for the chirality phase transition of Ising univer- 
sality. The natures of various possible Ising excitations 
are clarified in the dual model. 
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FIG. 1: Representation of the original triangular lattice (full 
line) and the dual hexagonal lattice (dotted line) sites labeled 
by lower case and upper case letters, respectively. Nearest- 
neighbor vectors are denoted as e± (Ei) through ez (E3) for 
the triangular (hexagonal) lattice. 




FIG. 2: The thin lines forming the triangular lattice have 
the chirality cr = 3 for upward triangles and cj = — 3 for 
downward triangles. The thick lines have the local chirality tij 
reversed from its ground state value, forming "defect bonds." 
The lines (green dotted) that intersect such defect bonds can 
form a closed loop (the smallest such loop is a hexagon) or an 
open segment. Each point of the line segment carries vorticity 
changes of ±<5 if it is connected on both sides and ±(1/2 — 5/2) 
if it is connected on one side only. 



